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Presented here are tables of values of % n and \„ , quantities defined by the 
eigenvalue problems 

(1 - x 2 )W ~ 2a*/ + (Xn - cV)*. = 

and 

\*.W - /' shl f x ~ f *Mdy. 

J -i r{x - y) 

In addition, some approximations for these quantities are given and evalu- 
ated. 



The prolate spheroidal wave functions of zero order, rp n (x), 
n = 0, 1, • • • , are bounded continuous solutions of both the differential 
equation 

and the integral equation 

KUx) = /' sin , c(x ~ » } Uv)dy. 

J-i t(x - y) 

The importance of these functions and the corresponding eigenvalues 
Xn and X„ for a great variety of problems, dealing with such diverse 
matters as lasers, communication theory, optics, noise theory, etc., can 
be found in the bibliographies of Refs. 1 and 2. It is our purpose here in 
response to numerous requests to present some numerical values for these 
eigenvalues. 
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Tables I and II list values of x« and X n respectively for n = (1 )20 (5) 
40 and c = (1 )20 (5 )40. The values given are, we believe, accurate to all 
eight figures listed.* Results of the computation are shown graphically 
on Figs. 1-4. 

The values of % n were obtained using the method of Bouwkamp as 
explained for example in Flammer. 1 This computation also gives expan- 
sion coefficients d r on (c ) in Flammer's notation, from which his quantity 
Z2 „ (1) (c,l) can be computed. The X's where then found from 

x« = -[fl „ (1 W)] 2 . 

T 

The tables presented required 0.027 hours of computing time on the 
IBM 7090. 

The following formulae for A„ and Xn are given in Ref. 2. For fixed n 
and small c 



-i[ 



2""(nl)' T , n+1 



(2n)!(2n + 1)1_ 
For fixed n and large c 

Some values computed from the terms explicitly exhibited in (1) and 
(2) are shown as dotted lines on Figs. 3 and 4. 

For n and c both large, we have the following result. Let b be fixed and 
let 



= r?(c + 61n2V3)] 



(3) 



where the brackets denote "integer part of". Then 

lim\„ = (1 + e ri ) -1 . (4) 



* The notation E ± XY following an entry in the tables indicates that the 
entry is to be multiplied by 10 ±xr where XY is an integer in decimal notation, e.g., 
E + 03 denotes a factor of 10 3 . 
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Fig. 1 — Eigenvalues, Xn , of (1 - z 2 )tf<" - 2a^' + (x - c 2 x 2 )^ = 0. 

The derivation of (3) and (4) given in Ref. 2 suggests the approximate 
formula 



X„ « X„ = (1 + e* 6 ')" 



6 = 



7T . 7T 

W 2" C + 4 



(5) 



(6) 



( 7 /2) +2 In 2 + £lnc 

for the near vertical rise portions of the X„ curves shown on Fig. 2. 
Here, y = 0.5772156649 • • • is the Euler-Mascheroni constant. The 
remarkable accuracy of this approximation is shown on Fig. 5. Here, for 
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Vn 0.5 




20 22 24 26 



Fig. 2 — Eigenvalues, X„ , of integral equation. Linear scale. 



fixed values of n and 6, we have determined values of c from (6) and 
for these values of n and c have plotted | (X„/X„) — 1 | vs n. It is seen 
that for 0.2 ^ X n ^ 0.9, (5) and (6) give an excellent approximation 
even for small values of n. 

Corresponding formulae for the Xn follow. For fixed n and small c 



ID 



x " = ni "' + 1} f 2 ! ' + (2n - \)(2n + 3) 



] c 2 + 0(c 4 ) 



and for fixed n and large c 



Xn = (2n + l)c - 



2n 2 + 2w + 3 



(2n + l)(w 2 + n + 3) 



16c 



•©■ 



+ -, 



If n and c become large according to (3) with b fixed, 



2 . n , .6—1 6—6 
x„ = c + The + -j g^ + 



»©■ 
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Fig. 3 — Eigenvalues, X„ , of integral equation for c < mr/2. Dashed lines are 
approximation (1). 






1758 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1965 



s 

6 
4 


- 


' \ VI \ 


\ \ \ V \ 




\ 




- 


V 


\ i 


\ 


_\_ 


\ 






\ 




- 


\\ 






\ 


\ 








\ 




_ 








\ 


\ 






I \ 




2 
10-2 

s 

6 

4 

2 

10-3 

a 

6 
4 

2 


n 


4_ 


i\ 


2 \ 3 \ a 


\ A 


\ 6\ 


7\ 


8\ 


A ,o \ 






\t 








\ 


\ \ 
\ I 










- 


\\ 








\ \ 










- 


y 










\ \ 
\ \ 




















\ \ 














ii 




\ 
























\' 




\ 
» \ 
\ \ 
\ \ 
















\ 


\ 




>\ \ 
















\ 


\ 




,\ 
















\ 


\ 




\\ 


















\ 




H 




\ ^ 












\ 


\ 








\ 




I0"« 

8 
6 

4 

a 












' 




\\ 




\ 












\ 






\l 




\ 












t: 




i 


V 




\ 










, 


\ 




\ 


\ ] 


\ 
I \ 


\ 




- 
















,\ 1 
\ 


i\ 










\ 


\ 




\ 


\ 


\\ 






10-5 

6 

4 

2 








\ 


\ 




\ 


\ 


\\ 






- 






\ 


\ 




\ 


\ 


\» 






- 






zt 


ZI 




zt 


zt 


\v 






- 












\ 


\ 














\ 




1 






















\ 




\ 








m-« 










1 


\ 




\ 


\ 


V 1 





8 10 12 14 16 18 20 22 

C 



Fig. 4 — Eigenvalues, X„ , of integral equation for c > mr/2. Dashed lines are 
approximation (2). 
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Fig. 5 — Accuracy of the approximation (5) — (G) for the eigenvalues X„ . 
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